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On v^/c^ expansion of the Dirac equation with external potentials 
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The «Vc^ 

expansion of the Dirac equation with external potentials is reexamined. A complete, 
gauge invariant form of the expansion to order (1/c)^ is established which contains two additional 
terms, as compared to various versions existing in the literature. It is shown that the additional 
terms describe relativistic decrease of the electron spin magnetic moment with increasing electron 
energy. 
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I. INTRODUCTION 

A semirelativistic expansion of the Dirac equation ( up 
to v"^ /c^ terms ) is treated in almost all sources on quan- 
tum mechanics which include elements of the relativistic 
quantum theory. This expansion is of importance for 
problems involving electric and magnetic potentials for 
which exact solutions of the Dirac equation do not ex- 
ist. The important examples include electron behavior 
in atoms, molecules and solids in the presence of a mag- 
netic field. Surprisingly, various versions of the w^/c^ ex- 
pansion quoted in the literature vary strongly from one 
source to the other 20j. Almost all final forms are 
incomplete, quite a few are not gauge invariant. 

The purpose of this contribution is to establish the 
complete and gauge invariant form of the j (? expan- 
sion of the Dirac equation with electric and magnetic 
potentials and to interpret its physical content. 

There exist two ways to carry the /c^ expansion. 
One is to consider large and small components of the 
wave function for electrons having positive ( or negative 
) energies and to find the large component to the desired 
order of 1/c by iteration. In this method one has to take 
into account a changed normalization condition for the 
large components ( to the same order ) . The other way is 
to use the Foldy-Wouthuysen unitary transformation in 
order to eliminate odd operators in the Dirac equation to 
the desired order. In this method the normalization con- 
dition is taken into account automatically since a unitary 
transformation does not affect the normalization. Both 
methods, if carried out consistently, lead to the same re- 
sult up to the order (1/c)^ ( see de Vries )• In our 
procedure we use the unitary transformation as it is more 
systematic. 



II. RESULTS AND DISCUSSION 



The Dirac equation for an electron reads 

H = I3mc^ + U + , 



(1) 



where m is the rest electron mass, U = —eV is the po- 
tential energy and O = c a ■ Tl is the kinetic energy. In 



order to avoid misunderstandings concerning the order of 
1/c, we put explicitly c in Eq. (1) and write the canonical 
momentum 11 = p+eA, i.e. the vector potential appears 
without c ( SI units). Symbol e signifies absolute value 
of the electron charge. Matrices /3 and a are taken in the 
well known Dirac form. We assume that the potentials 
V and A do not depend on time. 

The kinetic term O represents an odd part of the 
Hamiltonian mixing upper and lower components of the 
wave function. In order to eliminate the odd part the fol- 
lowing Foldy-Wouthuyscn transformation is introduced 



H' = exp{iS)Hexp{—iS) 



(2) 



where S = {—i/2mc^)f30 .The exponentials in Eq. (2) 
are expanded into power series to the desired order of 1/c, 
taking into account that S ^ 1/c. The choice of S elim- 
inates the odd term O in Eq. (1), but introduces other 
odd terms of higher order in 1/c. These are eliminated 
by two consecutive transformations S' and S" leading to 
the Hamiltonian , which is free of odd operators to or- 
[1/c^) . We do not quote this well known procedure 
The final result is 
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(3) 

where the brackets symbolize commutators. 

Evaluating the above quantities we will use a 4x4 spin 
vector operator 0, 0| 



S = 



cr 
a- 



(4) 



It has the following property for any two vectors C and 
D 

{a-C){a-D) = C-D + ii:-{CxD) (5) 
Using this identity we obtain 

~ mS-B , (6) 
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7^ oO^ = 

2mc^ 2m 



where /i = eh/ 2m is the Bohr magneton and B = V x A 
is a magnetic field. Clearly, the operator in Eq. (3) 
is just the square of given above. 
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The last term in Eq. (3) is calculated in two steps 

[O, C/] = [c a • n, [/] = chei{a ■ E) , (7) 

where E = — is an electric field. To evaluate the 
final commutator in Eq. (3) one needs two properties of 
the spin vector operator S, which follow directly from 
the corresponding properties of cr. We have 



OtkOij 



Skj + iskji^i , 



(8) 



(9) 



where dkj is the Kronecker delta and Sjki is the anti- 
symmetric unit tensor. The summation convention over 
repeated subscripts is employed. Using Eqs (8) and (9) 
we calculate 

[a n,a E]= aka-j [Ilfc, + aj]Ejnk 



/ c • ^ ^ ^dEj 



2iejki'>^iEjn 



= -ih{V ■E) + h^-{V xE)- 2iT, -{ExU) . (10) 

If the scalar potential V{y) has continuous first deriva- 
tives, then 'V X E = —V x 'VV — 0, and the second term 
vanishes. 

There exists an alternative way to evaluate the above 
commutator using directly Eq. (5) 

[a-n,a-E] ^ -ih{V-E)+ii:-{nxE-Exn) . (11) 

Both above forms are clearly equivalent. In agreement 
with the common practice we use the gradient sign V to 
emphasize that the differentiation concerns the electric 
field alone and not the wave function. 

The Hamiltonian (3) factorizes into two 2x2 blocks for 
upper and lower components of the wave function. For 
the upper block ( positive energies ) /3 is to be replaced 
by -1-1 and S by cr [cf. Eq. (4)]. This finally gives for an 
electron with positive energies the following Hamiltonian 
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mc 



eV- 



1 

2^' 



ficr ■ B- 



2mc^ 2m 



Am?c- 



:(T-{ExIi) 



:{V-E) 
(12) 

This is the main result of our paper. Now we briefly 
discuss its physical content. 

The fourth Pauli term is nonrelativistic although it 
results from the relativistic Dirac equation. As shown by 
Huang 1^ and Feshbach and Villars [l^], this term is 
related to electron's Zitterbewegung. 



The sixth term represents the spin-orbit interaction. It 
is written in the form calculated in Eq. (10) accounting 
for V X = 0. The alternative form is [ see Eq. (11)] 



eh 



:(T ■ {E xU-Ux E) 



(13) 



In both forms the canonical momentum 11 appears, as- 
suring the gauge invariance of the Hamiltonian (12). Un- 
fortunately, many well known references give noninvari- 
ant form with p. As demonstrated by Feshbach and Vil- 
lars [22], the spin-orbit term results from the linear con- 
tribution Ar to the electron displacement caused by the 
Zitterbewegung. 

The seventh term proportional to V • £J is the Dar- 
win term. It is commonly interpreted as coming from 
the quadratic contribution (Vr)^ to the electron displace- 
ment caused by the Zitterbewegung. 

Finally, we want to discuss the fifth term in the Hamil- 
tonian (12). Of all the enumerated references [1-20], only 
Corinaldesi and Strocchi , Messiah [l3] and Hecht po| 
quote the term (cr • II)'*, but, since it is nowhere sepa- 
rated into the orbital and spin parts and interpreted, it 
merits attention. To facilitate the discussion, it is help- 
ful to consider first a simple situation of an electron in a 
constant magnetic field B. In this case the Dirac equa- 
tion has exact solutions and the positive eigenenergies 
are MM 



[{mc^f + 2mc^D{n,p^,±)Y/'^ , (14) 



where 



D{n,p,, ±) = huj,{n +h) + 2l.±^B , (15) 

in which n — 0, 1, 2, ... is the Landau quantum number 
and LJr. = eB/m is the cyclotron frequency. Expanding 



the square root for 2D <^ mc^, one obtains 



D 



1 



2mc^ 



(16) 



We can identify the above expression with the first, third, 
fourth and fifth terms of the Hamiltonian (12) because, 
for the free electron, the eigenvalue of the orbital term 
U^/2m is hiJcin + 1/2) + pi /2m. Thus the fifth term 
in Eq. (12) corresponds to (1/c)^ order in expansion of 
relativistic energy (including the Pauli term /i cr ■ B). 

If we keep = but have otherwise arbitrary time 
independent magnetic field B(r), one can transform the 
Dirac Hamiltonian for positive energies into the form ( 
see Case Eriksen and Kolsrud 



2 n2 

H' = Pmc'[l + ^i^ 
mc Zm 



fi cr 



B)]l/2 



(17) 



Expanding the square root and retaining the first three 
terms we obtain exactly the corresponding terms in Eq. 
(12). 

Equation (15) shows that the spin contribution to the 
fifth term is not negligible in comparison to the orbital 
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contribution. In fact, for the n = Landau level the 
two contributions are exactly equal. In other words, the 
spin splitting of electron energies is equal to the orbital 
splitting. 

Turning to the physical meaning of the term in ques- 
tion, it is often stated that the — (l/8m'^c^)n* term in the 

j(? expansion reflects relativistic increase of the elec- 
tron mass with increasing energy. As to the spin term, it 
was demonstrated by Zawadzki that in the limit of 
vanishing magnetic fields the spin magnetic moment of a 
free relativistic electron is 



where TO(e) = e/c^ is the energy dependent relativistic 
mass. Thus, for the electron at rest, e = mc^, the spin 
magnetic moment reduces to the Bohr magneton, but, 
as the energy and the mass increase, the moment /i(e) 
decreases tending to zero. In the presence of magnetic 
field the decrease of the magnetic moment means that 
the spin splitting of the energy decreases with increasing 
energy. This decrease of the spin splitting can be seen 
from Eqs (14) and (15): for a given value of B both the 
orbital and the spin splittings diminish as n grows. It 
is then clear that the appearance of the complete fifth 
term in Eq. (12) which, as shown above, corresponds to 



the expansion (16) or, more precisely, to the expansion 
of Eq. (17), expresses not only relativistic increase of the 
electron mass but also relativistic decrease of the spin 
magnetic moment. 

Since the orbital and spin contributions to the fifth 
term commute, one can perform the indicated squaring 
directly. In this form the jt? expansion of the Dirac 
equation, as given by Eq. (12), contains two additional 
terms compared to the expressions given in the literature. 

In summary, we have critically examined the semirel- 
ativistic expansion of the Dirac equation with scalar and 
vector potentials. The complete, gauge invariant form of 
the expansion to order (1/c)^ is established. This form 
contains two additional terms, as compared to different 
expressions given in original papers and textbooks. It 
is demonstrated that the additional terms describe rel- 
ativistic decrease of the electron spin magnetic moment 
with increasing electron energy. 
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